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Abstract. The role of gradient dependent constitutive spaces is investigated on the exam- 
ple of Extended Thermodynamics of rigid heat conductors. Different levels of nonlocality are 
developed and the different versions of extended thermodynamics are classified. The local 
form of the entropy density plays a crucial role in the investigations. The entropy inequality 
is solved under suitable constitutive assumptions. Balance form of evolution equations is 
obtained in special cases. Closure relations are derived on a phenomenological level. 
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1. Introduction 

Weakly nonlocal thermodynamic theories are those that introduce the space derivatives 
of the basic variables into constitutive functions [Q. Second Law restricts considerably the 
form of the constitutive quantities and gives a genuine insight into the structure of the the- 
ories. Weakly nonlocal constitutive functions are mostly investigated in relation of material 
microstructure in mechanics (21 E] or to find nonlocal extensions of classical theories pQ. 

In this paper we investigate nonlocal constitutive spaces with different levels of nonlocality, 
namely of different order of derivatives. However, we derive also the restrictions that are due 
to locality assumptions on different levels. In our analysis we assume a non-equilibrium 
entropy function that can be approximated by its values measured at the equilibrium. Such 
an assumption is referred to as local state hypothesis [I]. 

We restrict ourselves to extended thermodynamic theories of rigid heat conductors El E| 
and introduce the heat flux together with a second order tensor as internal variables. The 
balance structure of the evolution equations is not postulated. Furthermore, the entropy 
current is regarded as a constitutive quantity and we are to give a complete solution of the 
thermodynamic constraints i.e. both the equalities and the residual dissipation inequality. 

In a previous work |Sj the local theory has been developed in the details and the evolution 
equations for fluxes of higher tensorial order have been obtained. Also it was proved that 
under particular assumptions on the entropy density and the entropy current the balance 
form can be recovered. Moreover, the system of equations was closed, in that the evolution 
equations for the highest order variable in the hierarchy - ordinary differential equations - 
can be interpreted as a closure relation. 

In the present paper we extend our investigation to the case of weakly nonlocal constitutive 
state spaces. The solutions are derived with the help of minimal assumptions on the form 
of the entropy, on its flux or on the evolution equations of the internal variables. In this 
way the different solutions of the entropy inequality are clearly classified. The most general 
assumption, that covers all existing phenomenological theories lies on the concept of current 
multipliers, which represent some constitutive functions entering the entropy current. We 
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will see that the final evolution equations are more general than the traditional balance ones. 
The conditions to recover the classical cases are clarified. 

In Section 2 we investigate first order nonlocality by applying Liu procedure [HI E3 f° r the 
exploitation of Second Law. 

In Section 3 we solve the Liu equation in the case of local state and local evolution equations 
for the dynamic variables. These assumptions result in a set of rather unusual restrictions 
from which we argue that some kind of nonlocality, either of the constitutive space or of the 
evolution equations, seems to be unavoidable. On the other hand, if we face with nonlocal 
state but local evolution equations then the entropy current is local, provided that the entropy 
density does not depend on the gradient of the internal energy. 

In Section 4 we investigate the traditional assumptions of Extended Irreversible Thermo- 
dynamics based on the following form of the local entropy jllj 

(1.1) s(e, q i} <$>ij) = s - ^mijqiqj - -Uiju^ij^kh 

where the matrices rriij and riijki are constitutive functions and sq is the equilibrium entropy, 
that depends only on the internal energy. We first suppose the entropy current is given as 

m 

(1.2) ji = Aijqj + B ijk $ jk , 

where A^ and Bij k are constitutive functions, the so called current multipliers. Then we 
explore the less general case ^Hl, too 

ds ds 

(L3) * = + 

It is worth noticing that the general form of the entropy current (jl.2|) reduces to (|1.3j) when 
Aij = §§<% and = J^r<%. We investigate different assumptions that can be compatible 
with the balance form of the evolution equations. 

In Section 5 we consider second order nonlocality but conserve the form (jl.lj) of the entropy 
density and the expression (|1.2j) of the entropy current. In such a case, due to the enlargement 
of the state space, the balance form can be obtained even if the general constitutive equation 
(|1.2|) holds true. We show that all previous examples can be recovered under simple special 
assumptions. 

In Section 6 we point out some nonlocal effects arising in thermal wave propagation at 
low temperature, which are described by the celebrated Guyer-Krumhansl equation ^2J ITB} 
116) . Such an equation has been derived by the authors by solving a linearized Boltzmann 
equation for phonon gas hydrodynamic. Here we prove that it can be obtained in the classical 
macroscopic framework of nonlocal irreversible thermodynamics. 

The previous results are discussed in Section 7, where a table shows the connections be- 
tween the constitutive assumptions and the thermodynamic restrictions, as far as the locality 
and nonlocality are concerned. 

2. First order nonlocality - exploitation of the Second Law 

In a rigid heat conductor at rest we start from the following local balance equation of the 
internal energy 

(2.1) e + q Vi = 0, 

where e is the density of internal energy, q\ i = 1,2,3 are the components of the heat flux, 
/ = //j = J||-, Xi i = 1,2, 3 are the Cartesian coordinates of the points of the body and 
the Einstein convention of summation over the repeated indices has been applied. The only 
equilibrium variable will be the internal energy e, while the first dynamic variable is supposed 
to be the heat flux c/j. As a further dynamic variable let us choose a second order tensor, 
whose components will be denoted by i,j = 1,2,3. Hence the basic state space (the 



wanted fields) in our investigations will be spanned by the variables (e, qi, This is a 13 
field theory, because the number of the independent fields is 13. However, some reductions 
are possible. For instance in Extended Thermodynamics tensor $jj is identified with the 
symmetric momentum flux Nuj-i coming from kinetic theory [7j. Then the unknown fields 
reduce to 10 and we face with a 10-field theory. Furthermore it is possible to decompose iV^j 
into an isotropic part, which is related to the internal energy, and a deviatoric part according 
to the equation jVjy] = \eb~[ij] + N < ij > , where N < ij > is symmetric and traceless ^7]. In such 
a particular case the internal energy coincides with one of the six independent components 
of Nun and we deal with a 9-field theory. We are investigating a first order weak nonlocality 
in all variables of the basic state, therefore the constitutive space is spanned by the basic 
state and its spacial derivatives, that is the fields (e, e/j, qy^ $ifk)- We assume that 

the evolution equations for the heat current qi and for $jj can be written in the following 
rather general form 

(2-2) en = g l7 

(2.3) <i>,. ; = f ij: 

where gi and fij are constitutive functions. With the assumption of first order nonlocality 
the spacial derivatives of the above equations give further restrictions |18[ I19j. 

(2.4) e H + q fji = 0, 

(2.5) q Vj - gij = 0, 

(2.6) ijj'fc - fif k = 0. 

These equations are sometimes referred to as prolonged forms of the evolution equations (|2.1|) , 
(tOll and (EU). 

The local balance of entropy is given by 

(2.7) s + ji'i = a s , 

with s standing for the entropy density, ji i = 1,2,3 for the components of the entropy 
current and o~ s for the density of entropy production. Second Law of Thermodynamics forces 
a s to be nonnegative. 

In the following we will investigate the restrictions from the inequality of the Second Law 
with the general assumption that both the entropy and the entropy flux are constitutive 
quantities. The method of the exploitation is given by the Liu procedure 9 : . However, 
according to our calculations, in the present case a generalized Coleman-Noll j2U] procedure 
gives identical results. 

Let us introduce the Lagrange-Farkas multipliers El F 1 , Tf and r|- for the evolution 
equations (|2.1I) . (|2,2j) and H2.3|) respectively. The multipliers Tf, Tfj and T^ k are for the 
prolonged evolution equations (j2.4|) . (|2.5|) and (|2.6j) respectively. 

Now, Liu procedure gives 

d 1 se + (d 2 s)iqi + (d 3 s)ij<&ij + (^s^e/j + (5 5 s)y-gj/j + 
+(d6s)ijk$ij'k + (diji)e>i + (d 2 ji)jqjH + (d 3 ji)jk^jk'i + {diji)j(i'ij + 
+{ddi)jkqj'ki + (dQji)jki&jk'u - 

I ' (e + ft'i) " I "7 (Qi ~ 9i) ~ r| (<i',, - fij) - Vj (<■', + qyp) - 
-r|- (qi'j - (dig^&j - [d 2 gi)kqk'j - (dsgiju^ki'j ~ {dm)ke*kj - 

— {d^gi)kiqk'lj — (d69i)klm^kl'lmj) 

~^ijk O&ij'k - (<hfij)&l - (9 2 fij)iqi'k - {dzfij)lm^lm'k - {d<ifij)l<i'lk - 

— (d5fij)lmQl'mk ~ {d6fij)lmn^lm'nk) > 0. 



Here d n , n = 1,2,3,4,5,6 denotes the partial derivatives of the constitutive functions 
according to the variables (e, qi, e>j, q^j, <&ij> k ) respectively. After some rearrangements 
of the inequality one obtains the Lagrange-Farkas multipliers from the first set of the Liu 
equations. These are obtained by imposing the coefficients of the time derivatives vanish. 

(2.8) dxs = T\ 

(2-9) {d 2 s) % = Tl 

(2.10) (d 3 % = r?-, 

(2.11) (d^ = If, 

(2.12) {d bS ) i3 = i*, 
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(2-13) (d 6 s) ijk = T b ijk . 

The second set of Liu equations is obtained by taking equal to zero the multipliers of the 
second order space derivatives. By applying (|2.8|) - (|2.13|) one can write them as 

(2.14) (d<iji)j + {d 5 s)ii(d 4 gi)j + (d 6 s) k ii(d A f k i)j = 0, 

(2.15) {dsji)jk - (d4s)i5 jk + (d 5 s) mi (d 5 g m )j k + (d 6 s) mn i(d 5 f mn )j k = 0, 

(2.16) (deji)jM + (d 5 s) mi (deg m ) jk i + (dQs) mni (d 6 f mn ) jk i = 0. 
Finally the residual dissipation inequality can be written in the following form 

[diji + {d 5 s)jidigj + {d 6 s) jk id 1 f jk ]e>i + 
+[(<hjj)i ~ d 1 s5 ij + (d 5 s) kj (d 2 g k )i + (d e s) k ij(d2f k i)i]qi'j + 
+ [{d?,jk)ij + {d$s) lk (dzgi)ij + {d 6 s)i mk (d 3 fi m ) ij }^ ij/k + 

(2.17) + (d 2 s) igi + > 0. 

It is easily seen that the Liu system (|2.14|) - (j2.16|) is composed by 117 differential equations 
constraining the set of the 832 partial derivatives of the constitutive functions s,ji,gi,fij 
with respect to the elements of the constitutive space (e, $y, e>i, q^j, &ij' k ). Without some 
simplifications there is no hope to solve such a system. 

In the following sections we are looking for special simplifying assumptions to solve the Liu 
equations (|2.14f) - (j2.16|) and the dissipation inequality (|2.17|) . First we will investigate cases 
where some of the constitutive functions are assumed to be local. 

3. Solutions of the Liu equations - locality assumptions 

3.1. Local state. Let us start the assumption of the local state in the form that the entropy 
is independent of the gradients: 

(3.1) s := s(e,qi,$ij) 

In this case the Liu equations ()2.14j) - (|2.16[) are simplified considerably and as a solution 
we obtain a local entropy current 

(3.2) ji= ji(e,qi,$ij). 
The dissipation inequality (|2.17|) simplifies, too 

diji&i + [{d 2 jj)i - disSi^qvj + {d 3 j k )ij^i jlk + (9 2 s)i9i + (d3s)ijfij = 

(3.3) = - disqin + (d 2 s)igi + (^s)^/^ > 0. 

3.2. Local state and local evolution. Let us investigate now the case when the evolution 
equations of the internal variables are ordinary differential equations, that is we assume that 
the constitutive quantities gi and fa take the local form: 

(3.4) gi := gi(e,qj,$ jk ) 
(3-5) fa := fij(e,qj,$ jk ) 



Now the dissipation inequality simplifies further the possible constitutive functions, because 
the coefficients of the derivatives should disappear. These restrictions result in a rather 
unusual material, since we get 

dj' ds 

(3.6) j, W . = g-Jir 

On the other hand the dissipation inequality can be written in a force-current form and 
can be solved for the constitutive functions gi and 

(d2s)igi + (d 3 s)ijfij > 0. 

By (|3.6|) it follows that the temperature of the material g| is independent of the internal 
energy. Such a property is in contrast with the physical reality. We conclude that some 
nonlocality is necessary in modelling rigid heat conductors. 

3.3. Local evolution. Let us assume now that the evolution of the internal variables is local, 
but there is no local state, therefore we require (|3,4j) and (|3.5|) . but (|3.1|) is not assumed. In 
this case the entropy current is nonlocal, but the nonlocality is rather reduced. The Liu 
equations (|2. 14|) and (|2.16|) give that the entropy current does not depend on the gradients 
of e and and (|2,15|) simplifies to 

dji ds 

The nonlocality in the % is due to balance form of the evolution equation of the internal 
energy (|2.1|) . From the above constraint one can easily see, that the entropy current is local 
if we further assume that the entropy is local in the internal energy, 

s := s(e,qi,$ij,q V j,$ij> k ). 



4. Solutions of the entropy inequality in case of local state 

As we have seen above, the Liu equations are trivially solvable in the local state. How- 
ever, the solution of the dissipation inequality can be achieved only with further assumptions. 
Moreover, there are different assumptions to have physical models, to introduce a suitable 
gradient dependencies. In Classical Irreversible Thermodynamics |5, and in Rational Ther- 
modynamics |2U| the mentioned requirement of nonlocality is achieved by introducing the 
gradient of temperature (or, equivalently, of the internal energy) into the constitutive space. 
In Extended Thermodynamics 00 the constitutive space is local but the evolution equations 
are balances, they have a special nonlocal form. All kind of theories of Extended Thermody- 
namics resulted in more of less satisfactory models of different phenomena, but they provide 
different solutions of the entropy inequality. 

4.1. Local state and special nonlocal evolution: linear nonlocality. In this case the 
evolution equations depend linearly on the gradients: 



(4.1) gi :— AijCj + Bij k q k /j + Cijki&ki'j, 

(4.2) fij := Dijke/k + Eijkiqi'k + Fijkim&im'k, 

where Aij,B ijk , C ijk i, D ijk , E ijkU FjjMm are local constitutive functions. 
Now the dissipation inequality 1)2. 17|) reduces to a solvable form as 

(4.3) (diji + (d 2 s)jA jk + (d 3 s) jk D jki )eii 

(4.4) +((<9 2 jj)j - disSij + (d 2 s) k B kji + (d^s) k iE k iji)q v j 

(4.5) +({ddk)ij + (d2s)iCi kij + (d 3 s)i m Fi mkij )^ ifk > 0. 



As the quantities in the parentheses are local functions they should be zero respectively. 
Therefore we get 

(4.6) diji = -{d 2 s)jA jk - (d 3 s) jk D jki 

(4.7) {d 2 jj)i = dtsSij - {d 2 s) k B kji - (d 3 s) k iE k i ji 

(4.8) {dzj k )ij = -(d 2 s)iCi kij - (d 3 s)i m Fi mkij . 

These equations cannot be solved without any further ado. However, we can see that even 
if we do no know anything on the entropy current ji they result in strong correlations on the 
entropy derivatives and the evolution equation, as the mixed partial derivatives of ji should 
be equal. On the other hand let us observe that in this case the entropy production is zero, 
there is no dissipation. 

4.2. Local state and special local evolution: balance form. In this case one assumes, 
that the evolution equations have a special balance form. Therefore there are potentials Qij 
and H ijk of the fields (Ay, B ijk ,C ijk i)(e,qi,®ij) and (D ijk , E ijkh F ijk i m )(e,qi,$ij) respec- 
tively. The evolution equations can be written as 

(4.9) Qi := diQij&j + (d 2 Qij) k q k ij + (dsQ^ki^kl'j, 

(4.10) fij := diHij k e/ k + (d 2 Hij k )iqii k + (d 3 Hij k )i m ^i k / k . 

The conditions (|4.6[) - (|4.8|) can be written as 

(4.11) diji = -(d 2 s)jdiQij - (d 3 s) jk d 1 H ijk 

(4.12) (d 2 jj)i = dis5ij - (d 2 s) k (d 2 Qij) k - (d 3 s) k i(d 2 H ijk )i 

(4.13) (d 3 jk)ij = -{d 2 s)i(d 3 Qij) k i - (d 3 s)i m (d 3 H ijk )i m . 

As a consequence the above system of equations can be solved, as the entropy current is a 
potential of the field (qi,Qij, Hij k ), therefore it can be conveniently written as 

(4.14) ji(e, q h $y) = ji(qi, Qij,H ijk ). 

This fact can be expressed also with differential forms, according to the traditions of 
thermodynamics 

(4.15) dji = disdqi + (d^jdQji + (d 3 s) jk dH jki = rM^ + TfdQ^ + T) k dH jki . 

The derivatives of the entropy current are identical to the intensives, the derivatives of 
the entropy functions. However, the variables are different. This form results in serious 
restrictions of the entropic intensives, and the currents Qij and Hij k because the mixed 
second partial derivatives should be equal in these variables, too. Expressed in the basic 
variables the above requirements are rather ugly 

(4.16) {d 22 s) kj diQ k i + {d 23 s)i kj d 1 Hi ki = d n s5ij + (d 12 s) k (d 2 Q ik )j + {di 3 s)i k (d 2 H m ) j , 

(4.17) (d 32 s)i jk diQu + (d 33 s)i mjk diHi mi = {di 2 s)i(d 3 Q u ) jk + (d 13 s)i m (d 3 Hi mi ) jk , 

{d 3 is) k i5ij + (d 32 s) mk i(d 2 Q im )j + (d 33 s) 

mnkl 

(d 2 H mn i)j — 

(4.18) = (d 22 s) m i(d 3 Q m i) jk + {d 32 s) mn i(d 3 H mni ) jk . 

The property (|4.15|) is an important consequence of the balance form of the evolution 
equations. It is independent of the choice of the basic variables. If one assumes e.g. that the 
chosen internal variable is the current of the heat flux, = Qij, as it is usual in extended 
thermodynamics, then the above system of requirements simplifies but does not disappear. 

In Rational Extended Thermodynamics it was shown that the above result of the phe- 
nomenological theory can be in accordance with the kinetic theory of gases, at least with a 
classical formulation of kinetic physics. A crucial step in the different systems was the choice 
of the phenomenological variables (we will see, that all of the currents cannot be chosen as 
internal variables without any further ado) and the use of source terms in the balances. 



4.3. Local state and balance form evolution: isotropy. The system (|4.16j) - (j4.18j) does 
not have a general solution for the currents. Qij and H^^ cannot be determined by the 
entropy function in general. Therefore we have lost one of the basic flavors of irreversible 
thermodynamics, that the requirements of the second law can be exploited constructively to 
find a the appropriate evolution equations. Now the dissipation inequality was solved, but 
the evolution equations cannot be determined constitutively. 

Jou, Lebon, Mongiovi and Peruzza gave some simplifying conditions to have a solution of 
the conditions (|4.16j) - ()4.18|) on the phenomenological level |21j . They have assumed a local 
state, balance form evolution equations and a simpler set of variables, they introduced only 
qi as an additional variable. Then only condition (|4.16j) applies in a simplified form as 

(d22s)kjdiQki = dnsSij + (di 2 s) k (d 2 Qik)j- 

Moreover they have assumed isotropic materials, when all scalar valued functions, including 
the entropy, depend only on q 2 = qiqi and the flux of the heat current and the entropy current 
can be written as 

(4.19) Qij = p(e,q 2 )+iP(e,q 2 )q iqj , 

(4.20) ji = ^{e,q 2 ) qi . 

Now the requirement (|4.15|) results in the following system of equations, as (4.2)-(4.4) in 

(4.21) d e * = 2d q 2s(d e (3 + d e 7Pq 2 ), 

(4.22) # = d e s + 2d q 2 sipq 2 , 

(4.23) d q 2$> = 2d q 2s(d q 2(3 + d q 2il)q 2 ) . 

Therefore the entropy current can be written as 

(4.24) j i ={d e s + 2d q 2S^q 2 )q l . 

After further calculations, considering also the requirement ()4.16j) . one can get explicit 
solutions for the functions (3 and ip together with some additional restriction on the form of 
the entropy function. 

5. Local state and special entropy current 

A different solution of the dissipation inequality can be given with the help of the entropy 
current. As in our previous work [B] we consider the local entropy of the form (jl.ljl . This 
general form is motivated by the requirement of the thermodynamic stability, or, equivalently 
the requirement of the concavity of the entropy function on the non-equilibrium part of the 
state space (spanned by qi and Therefore, and riijki are positive definite constitutive 
functions. Moreover, let us assume the entropy current takes the form (|1.2|) and let us 
introduce the convenient notations m, n, h, rh as follows 



(d 2 s)i = -mijqj - -^{d 2 m lj ) i qiq : j - -{d 2 n rjk i)i$rj$kl 

(5.1) = -rhijqj - (d 2 n jk )i$ jk , 

(d3s)ij = -^(d 3 mi k )ijqiqk - n ijk i$kl ~ -jipznikrs)ij^lk^rs 

(5.2) = -(d 3 m k )ijq k - h ijk i§ kl 



If my and are constant, then rhij = rriij, hij k i = n^ki, {d 3 fhk)ij = and ( d 2 hj k )j = . 

Let us emphasize again that the entropy, written in the Gyarmati form and the 

entropy current, written in (|1.2|) are only convenient notations as long as the corresponding 
inductivities and current multipliers are general constitutive functions. 



With the (|1.2j) form of the entropy current and using the notations of ()5.1|) and (|5.2|) the 
dissipation inequality (|3.3j) in local state can be written in the following form: 

[Aji'j - rhijgj - (d 3 m) kj if ik ]qi + [Ay - dis5ij}qj>i + 

(5.3) [Bkij'k - (d2n)kji9k — riijki fki)®ij + Bij k &j k >i > 0. 

Seemingly the system is a normal force-current system, because the coefficients of the ther- 
modynamic forces qi, qij, &ij k all contain undetermined constitutive quantities Aij,gi, fij, 
Bijk respectively. However, let us observe that in local state gi and fij are nonlocal but all 
other terms are local in the above inequality. This fact introduces degeneracy since, although 
the coefficients of the derivatives q^j and $ij> k cannot disappear, their possible couplings are 
rather reduced, e.g. Bij k is local, therefore cannot depend on its own force $ij' k - 

In this degenerate case the solution of the dissipation inequality is not straightforward. 

Fortunately we can avoid the treatment of degeneracy, e.g. by assuming that rriij and riij k i 
depend only on the internal energy and introducing the form Hl.Mj) of the entropy current with 
the assumptions Ay = d\s5ij and By k = ($2s)fc<%- In this particular case the dissipation 
inequality reduces to 



(5.4) ((<9is)'j - rriijgj - m ki ® jklj ) qi + {-mjiq Vi - n ijk ifki) $ij > 0. 
This is a force-current system, with the following Onsagerian solution 

(5.5) - rriijiqj + ^kj'k) + (dis)'i = Ljjqj + L\? k <S> jk , 

(5.6) -riijki&kl - rrijiqi'i = L^j k q k + L^ kk ^ kh 



The system above conserves the structure already obtained in [Hj in the case of local state. 
One should emphasize the central role of the invertibility of the matrices rriij and Tiij k i in 
order to obtain the balance form. Such a property is not trivial since there exist real materials 
for which it is not guaranteed. A classical example is given in where an electric circuit 
described by dynamic variables is considered. 

As it was observed in [S] the second equation is a phenomenological closure of the system 
(|2.1|) and (J5.5|) . It is remarkable that we have recovered the usual phenomenological structure 
of Extended Irreversible Thermodynamics keeping the entropy current of the form (|1.3jl . The 
same is not true for the constitutive equation (jl.2j) which, in the case of first order nonlocality, 
seems to be much too general. However, the compatibility of the entropy current of Verhas 
(jl.3j) and the requirements of the balance form (|4.15|) is valid only with the restriction that 
the conductivities depend only on the internal energy. The mentioned solution of Jou et. al. 
(|4.24j) clearly does not have the Verhas form and indicates the necessity of a more general 
treatment. 

In the next section we study the same problem in the presence of second order nonlocality. 
6. Second order nonlocality - solution in local state 

Now we will extend our investigations to consider second order nonlocalities. However, 
for the sake of simplicity we will investigate only the case of local entropy, i.e. systems 
in local state. The basic state is spanned by the variables (e,qi,&ij) as previously. How- 
ever, the constitutive space contains the second order space derivatives and is spanned by 
(e, qi, &ij,e>i,qi'j, &ij' k , &y , qi>j k , §ij> k i)- Therefore, in the exploitation of the entropy inequal- 
ity (|2,7|) we should consider as constraints the evolution equations (|2.1j) . ()2.2j) and ()2.3|) . their 
first prolongations ()2.4j) . ()2.5|) and ()2.6|) and also their second prolongations as follows 

(6.1) e>ij + qf jik = 0, 

(6.2) qi>j k — gij k = 0, 
(6-3) &ij'ki - fij'ki = 0. 

Our simplifying condition of local entropy can be written as 

s = s(e,qi,$ij). 



Let us introduce again the Lagrange -Farkas multipliers T 1 , rf, r?., Tf, r?. and T 6 ijk for 
the evolution equations (|2~T|) . (J22J), <|2~H|) and their prolonged forms J23J), (|2~5|) and (|2~fi|l 
respectively. The multipliers Tjj, r 8 jfc and r? jfc/ stand for the second prolongations (j6.1j) . 
H6.2[) and ()6.3|) respectively. The Liu procedure gives 

dise + {d 2 s)icti + (d^ij&ij + 
+(diji)z'i + {d 2 ji)jqj'i + {d^ji) jk ^jk'i + (34j;)j e 'r/ + (d53i)jk<lj'M + 

+ {dQji)jkl^jk'li + {dl3i)jk&'ijk + {d%ji)jkiqj'kli + (dgji)jklm^jk'lmi 

i (e + - r? (ft - - r| (\- - /y) - rf + ^^) - 

- r |fe'i - (di9i)e>j ~ {d29i)kqk'j ~ {d<i9i)kl$kl'j ~ {dm)ke>kj ~ 

— {d59i)klQk'lj — {d§9i)klm$kl'mj ~ {d-jgi)kie'klj — 

— {d%gi)klm<lk>lmj — (pV9i)klmn®kl'mnj\ 

~^%k&ij'k — (9lfij)e>k - (d2fij)iqi'k — (d3fij)lm®lm>k ~ {difi^l&lk ~ 

— {dbfij)lm<ll'mk ~ ($6 fij) lmn$ Ira' nk ~ (^7 fij)lm e 'lmk ~ 

— (dsfij) ImnQl'mnk - (dgfij) Imno^ lm! nok\ 

— rjj (eij + qk'kij) — T%k [<li'jk - 9i'jk) ~ ^%kl (^ij'kl ~ fij'klj > 0. 

Here d n , n = 1,2, 3, 4, 5, 6, 7, 8, 9 denotes the partial derivatives of the constitutive functions 
according to the variables of the constitutive space (e, ft, <&y, en, q^j, ^ij'ki^'ijiqvjki^ij'kl) 
respectively. The first set of Liu equations defines the Lagrange-Farkas multipliers as the 
derivatives of the entropy and gives that the last six multipliers are zero, due to the local 



state. 
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x ijk 


= o, 
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r 9 

1 ijkl 


= 0. 



Considering l|6.7|) - (|6.12|) the second set of Liu equations are also simple, 

(6.13) (d 7 ji) jk = 0, 

(6.14) (d 8 ji)jki = 0, 

(6.15) (d 9 ji) jk i m = 0. 

Therefore the entropy current presents only first order nonlocalities 

ji = 3i(e, ft, ®ij,e>i, qi>j,®ij>k)- 

Considering the above conditions the dissipation inequality can be written exactly in the 
same form as it was in case of first order nonlocal constitutive space ([3 .3)1 : 

(6.16) (ji)>i - disqvi + (3 2 s)ift + (d 3 s)ijfij > 0. 

On the other hand now the constitutive quantities are higher orderly nonlocal, the entropy 
current is a first order nonlocal quantity and ft and fij are second order nonlocal quantities. 
Let us introduce the same assumptions on the form of the entropy and of the entropy current 



as previously with assuming and (|1.2|) with the notations ()5.1|) and (|5.2|) but letting the 
current multipliers and B^ k to contain first order nonlocalities. The dissipation inequality 
can be written in the same form as above 

[Aji'j - rhijgj - {d 3 m) kji f ik ]qi + [Ay - d 1 s5 ij }q jli + 

(6-17) [Bkij'k - {fanjkjiQk - hijkifki]&ij + Bijk&jk'i > 0. 

However, in this case it is a non-degenerate force-current system, due to the extension of 
the constitutive state space. All additive terms contain unknown functions. Therefore one 
can give an Onsagerian solution as follows 

- rhijgj - {d 3 m) ijk f jk + Aji'j = 

(6.18) Ljjqj + Ljf k q j/k + Lj'f k ® jk + L ijkl $j k n, 

(6-19) Aij - disdij = L 2 j k q k + L 22 kl q kn + L 23 kl $ k i + L 2 jklm $ kVm , 

— {d2n) k jig k — hij k if k [ + B ki y k = 

(6-20) L^j k Qk + L ff kl q k >i + L 33 kl $ kl + L 3 f klm $ kVm , 

(6-21) Bij k = L^j kl qi + L^j klm qi/ m + L^ 3 klm ^i m + L^ 4 fcZmn $/ m / n . 

Here L 11 , L 12 , L 13 , L 14 , L 21 , L 22 , L 23 , L 24 , L 31 , L 32 , L 33 , L u , L 41 , L 42 , L 43 , L 44 are constitu- 
tive functions with the suitable definiteness restrictions. Moreover (|6.18j) and (|6.2Uf) are the 
candidates of balances under suitable restrictions. It is worth noticing, that the current 
multipliers A^ and B, L j k are given explicitly. Therefore, they can be easily eliminated from 
the above system substituting (|6.19|) into (|6.18|) and (|6.21|) into ()6.2Uj) . The resulting set of 
equations is closed and contains first and second order derivatives of the basic state. 

- rhijqj - (d 3 m) ijk & jk + 

+ [disSij + L 2 j k q k + L 22 kl q k >i + L 23 kl $ H + L 2 jklm $ k i> m ]>j = 

(6-22) = L}}q j+ L} 2 k q fk + Llf k ^ jk + L}^ jkn , 

-(d2n) kji q k - n ijkl <& kl + 

J fjkl ( ll + ^ijklmMl'm ~+~ ^ijklm^lrn ~t~ ^ijklmn^lm'nVk 
4jkQk + L 3 fkl<lk>l + L 3 jkl^kl + Lij klm $ k ii m , 

As one can see, the evolution equations are rather general. They are more general that 
those given with the help of the entropy currents On the other hand the entropy 

current (|4.24l) is also a special case of the general p.2|) . However, the compatibility to the 
current potential structure, expressed by (|4.15|) . cannot be expressed explicitly in general. 
(|1.2|) is definitely more general regarding the nonlocality, moreover, ()4.24|) proves, that can 
be compatible with the potential structure in special cases. However, it definitely restricts 
the functional form of the current multipliers A^ and By k . 

Moreover, the potential structure is unavoidable requiring locality and balance form evo- 
lution equations. If rriij and nij k i are constant a balance structure similar to that obtained 
in [S] is recovered. This fact represents the main effect of the enlargement of the state space 
since the balance structure is compatible with more general entropy fluxes. 

7. Nonlocal second sound: the Guyer-Krumhansl equation 

Thermal wave propagation, sometime referred to as second sound, is a low temperature 
phenomenon which can be observed, for instance, in dielectric crystals such as Sodium Flu- 
oride (NaF) and Bismuth (Bi) |231 124| ITT] . It requires an extension of the classical Fourier's 
theory in order to remove the paradox of infinite speed of propagation of thermal disturbances 
[25] . Phonon gas hydrodynamics EJJ supplies a satisfactory explanation of heat transport 
at low temperature. Phonons are quasi-particles which obey the Bose-Einstein statistics. In 
a solid crystal they form a rarefied gas whose kinetic equation can be derived similarly to that 



+ [LijklQl + ^ijklmQl'm + Lij k i m <&lm + ^ijklmn^l 

(6-23) = L 3 } k q k + L 32 q kn +L 33 ^ kl +m klm ^ k 



of an ordinary gas. Phonons may interact among themselves and with lattice imperfections 
through two different types of processes: 

i) Normal-(N) processes, that conserve the phonon momentum; 

ii) Resistive-(R) processes, in which the phonon momentum is not conserved. 

The frequencies vn and vr of normal and resistive processes determine the characteristic 
relaxation times tn = and tr = Diffusive processes take over when there are many 
more R-processes than N-processes. If instead there are only few R-processes and many more 
N-processes, then a wave like energy transport may occur. 

Such a phenomenology is satisfactorily described by the Guyer-Krumhansl equation [Tl|, eq. 
(57) (see also 

(7.1) 4i + —qi + \(?c v T>i = \t n c 2 (qvjj + 2qj ,jA , 

tr 3 5 v 

where T = T — To is the temperature variation (T is the temperature and To is the average 
temperature), cy is the specific heat and c means the Debye phonons velocity. Such an 
equation, which generalizes the Maxwell-Cattaneo-Vernotte equation [23] 

(7.2) T R q\ + qi = -T R ^c 2 c v T>i, 

was the first in the literature to include both relaxation times. It can be obtained by the 
linearized Boltzman equation for phonons in the Debye approximation, if one maintains 
terms 0(tn) |T4"] . The material coefficients tn, tr and cy are all depend on the temperature. 
According to experimental observations JH] one can get 

cy 



aT 3 , T R = d e" a/T , t n = bT' 



Here a, b, d, a are constant coefficients and m € {3, 4, 5} depending on the material. Let 
us remark that the last function, a temperature dependent t^t is clearly contradicts to the 
assumptions made by Guyer and Krumhansl ^1] p771, and did not consider, that (|7.1|) is a 
linearized equation. 

In Extended Thermodynamics the Guyer-Krumhansl equation can be recovered in a 4-field 
theory provided one assumes weakly nonlocal constitutive equations for the internal energy 
e and for the momentum flux 17 . 

However, such an approach seems to be questionable because: 

i) a nonlocal constitutive space is in contrast with the basic assumptions of Extended Ther- 
modynamics. 

ii) a nonlocal internal energy does not assure that the specific heat cy = is positive. 

iii) a nonlocal internal energy would modify also the energy balance. 

The last observation was pointed also by Dreyer and Struchtrup and they suggested to 
consider higher order moments El- The other important question is the temperature 
dependency of the coefficients. A phenomenological theory cannot predict the exact form to 
the constitutive functions, but gives restrictions and interrelations. These restrictions on the 
temperature dependencies are frequently treated rather loosely to get the compatibility with 
the linearized kinetic theory [141 IB]. 

These problems do not arise in the present theory. In fact, we can obtain the equation (|7.1|) 
by considering a 4-field model with second order nonlocality, namely the balance equation 
(|6.22|) with &ij = and rriij depending only on the internal energy. It yields 



(7.3) -mijdj + [disSij + L ijk q k + L ijia qm\ = L ij Qj + L ijkQfk- 



It is convenient to choose as equilibrium thermodynamic variable the absolute temperature T a 

-21 _ rl2 
J ijk ijk 

1 _ ds 



instead of the internal energy e. Let us assume the constitutive equations Lff, = L 12 fc = 0. we 



may exploit the thermodynamic relation = §f and introduce the temperature perturbation 



T by T a = To + T, where To is a background (average) temperature. Now, equation (|7.3jl 
reduces to 

(7-4) mijqj + Lf/gj + ^ T/; = L^ kl q kqj . 

Here m, L 11 and L 22 can depend on the background temperature To. Finally, equation Q7.1|) 
is easily obtained by Q7.4J1 under the further constitutive assumptions 

3 

- 2 ^T 5 ' 



(7.5) rriij — „ 5 



(7-7) L?? w = —^L—SijSu. 

That way, the Guyer-Krumhansl theory may be obtained in a classical macroscopic frame- 
work. Let us remark, that the previous phenomenological linearization is only a last step in 
the calculation based on the linearization of the Boltzmann equation (eq. (57) to eq. (59) in 
|14j). A similar result has been obtained in |29j . 

8. Conclusions 

We have shown that Classical Irreversible Thermodynamics supplemented with dynamic 
degrees of freedom is consistent with the idea of higher order fluxes as independent thermo- 
dynamic variables. In such a framework the balance form of the evolution equations can be 
obtained under suitable constitutive assumptions. The locality of the entropy density, i.e. 
the local state assumption, plays a central role. Moreover, the nonlocality of the constitutive 
functions results in a wide class of materials including the classical Cattaneo's and Guyer- 
Krumhansl heat conductors, which normally are derived from kinetic theory. Let us observe 
that the model above encompasses all extended thermodynamic models, since the balance 
structure represents a very particular form of nonlocal evolution equations (|2.2|) and ()2.3|) . 

The figure below shows the connections between the constitutive assumptions and the 
thermodynamic restrictions, as far as the locality and nonlocality are concerned. The last 
column refers to the corresponding section of the paper. La denotes a locality assumption. 
The double head arrows denote specific interrelations, e.g. the black double head arrows in 
the first row denote equations ()2.14)) - (|2.16|) . that give conditions between the entropy and 
entropy current functions and evolution equations. The fifth row refers to Rational Extended 
Thermodynamics and the sixth row to the specific solution of the potential restriction given 
by Jou, Lebon, Mongiovi and Peruzza. As in Extended Thermodynamics the dissipations 
inequality is fulfilled as an equality, we cannot conclude anything on the dissipative consti- 
tutive functions (that is why there are no white arrows in these rows). Let us observe, that 
assumptions on the form of the entropy current makes possible to build all requirements of 
the Second Law into the evolution equations in general (last two rows). In this case the 
Second Law become a material property, satisfied independently of the initial conditions and 
the resulted constitutive functions depend only on the material. 

It is worth noticing that the potential form of the entropy current and the balance structure 
is compatible with the general form of the entropy current, such as that proposed by Nyfri, 
which yields (|1.3|) when = d\s8ij and Bijk = (c^s)/^-. In this case the balance structure 
results in restrictions for the current multipliers from the potential requirement (|4.15|) . 

We obtained closure relations both with local and nonlocal entropy current and in the first 
case the closure for the highest order flux was an ordinary differential equation. The obtained 
thermodynamic closure of the hierarchical structure stresses some deeper relations between 
the thermodynamic and the more detailed kinetic structure, similar that was recognized e.g. 
in [HI. 



We reinspected the thermal wave propagation at low temperature and proved that the 
Guyer-Krumhansl theory of second sound can be derived in the framework of macroscopic 
nonlocal irreversible thermodynamics. 
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Figure 1. Notations: L - local, nL - nonlocal, L(e) - local in the variable e, 
La - assumed locality, L, - i-th order nonlocal, Black arrow - consequence of 
the Liu equations, White arrow - consequence of the dissipation inequality 
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